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Abstract 

A  new  model  of  equilibrium  current  sheets  in  a  collisionless  plasma  incorporating  ion 
flows  that  are  asymmetric  and  sheared  across  the  current  sheet  is  developed.  Ions  are 
treated  as  single  particles  and  electrons  as  a  massless  fluid.  The  resulting  current  sheet  is 
a  Vlasov  equilibrium  satisfying  Ampere’s  law.  The  current  sheet  thickness  is  shown  to  be 
of  the  order  of  the  ion  Larmor  radius.  It  is  found  that  the  structure  of  the  current  sheet 
depends  on  the  distribution  function  of  the  ions  entering  the  current  sheet.  A  characteristic 
feature  of  this  class  of  equilibria  is  that  the  pressure  tensor  is  anisotropic  and  nondiagonal 
while  the  reversing  component  of  the  magnetic  field  is  similar  to  the  hyperbolic  tangent 
Harris  field.  The  results  show  that  the  asymmetric  sheared  ion  momentum  flows  as  well  as 
the  off-diagonal  pressure  tensor  elements  play  an  essential  role  in  the  force  balance.  The 
model  is  applied  to  current  sheets  eo-moving  with  the  solar  wind  and  is  directly  compared 
with  high-resolution  WIND  MFI  magnetic  field  data  and  WIND  3DP  particle  data.  The 
agreement  is  good.  The  theoretical  results  are  scalable  to  collisionless  plasmas  in  different 
regimes. 
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1.  Introduction 

Current  sheets  occur  over  a  wide  range  of  regimes  in  solar-system  and  astrophysical 
plasmas  including  the  magnetosphere,  the  solar  wind,  and  the  solar  corona.  Such  current 
sheets  can  store  a  significant  amount  of  magnetic  energy,  and  the  release  of  this  energy  can 
have  important  consequences.  For  example,  the  disruption  of  the  night-side  magnetotail 
current  sheet  is  thought  to  cause  geomagnetic  substorms  [e.g.,  Dungey ,  1961;  Russell  and 
McPherron,  1973;  Fairfield ,  1984;  Hones  et  al,  1984].  In  the  solar  and  stellar  coronae, 
current  sheets  are  assumed  to  play  an  important  role  in  coronal  heating  [ Parker ,  1984]  and 
solar/stellar  flares  [ Priest ,  1982].  In  laboratory  reconnection  phenomena,  current  sheets 
in  confined  plasmas  have  been  observed  [e.g.,  Yamada  et  al.,  2000;  Carter  et  al,  2002; 
Fumo  et  al.,  2006;  Lawrence  and  Gekelman,  2009]. 

Perhaps  the  best  understood  natural  current  sheet  is  the  magnetotail  current  sheet, 
which  has  been  extensively  studied  since  its  discovery  [Wess,  1965].  The  current  sheet  is 
observed  to  become  “thin”  in  association  with  the  occurrence  of  substorms,  characterized 
by  thickness  of  the  order  of  (several  times)  the  thermal  ion  Larmor  radius  based  on  the 
field  just  outside  the  current  sheet  [e.g.,  Fairfield,  1984;  Frank  et  al,  1984;  McPherron  et 
al,  1987;  McComas  et  al,  1986;  Mitchell  et  al,  1990;  Lui  et  al,  1992;  Sergeev  et  al,  1993; 
Asano  et  al,  2003;  Runov  et  al,  2006].  The  magnetotail  current  sheet  is  embedded  in  the 
thicker  central  plasma  sheet.  In  the  solar  wind  (SW),  current  sheets  are  encountered  in 
the  form  of  rotational  and  tangential  discontinuities  [e.g.,  Lepping  and  Behannon,  1986], 
Their  thicknesses  are  also  of  the  order  of  the  ion  Larmor  radius  [e.g.,  Siscoe  et  al,  1968; 
Burlaga  et  al,  1977].  The  heliospheric  current  sheet  embedded  in  a  sector  boundary 
may  be  one  order  of  magnitude  thicker  than  such  discontinuities  [e.g.,  Winterhalter  et  al, 
1994].  For  the  solar  corona  and  astrophysical  systems,  however,  no  in  situ  measurements 
of  current  sheets  are  available. 

Current  sheets  are  formed  in  boundary  layers  between  neighboring  regions  with  oppo¬ 
sitely  directed  magnetic  fields.  If  the  regions  have  nearly  the  same  properties,  the  current 
sheets  are  nearly  symmetric.  The  magnetotail  current  sheet  and  the  heliospheric  current 
sheet  are  two  such  examples.  However,  current  sheets  in  regions  such  as  the  magnetopause, 
the  solar  corona,  and  the  solar  wind  may  occur  between  regions  with  dissimilar  properties. 
They  are  expected  to  be  asymmetric.  Current  sheets  have  been  extensively  studied  using 
symmetric  geometries,  but  only  limited  theoretical  understanding  of  asymmetric  current 
sheets  exists  at  this  time.  In  the  present  paper,  we  present  results  of  a  study  in  which  a 
self-consistent  asymmetric  current  sheet  equilibrium  is  constructed  using  a  hybrid  (kinetic 
ion  and  massless  fluid  electron)  simulation  approach,  taking  into  account  the  exact  ion 
orbits  and  such  properties  as  the  anisotropic  and  nondiagonal  pressure  tensor  and  sheared 
ion  flows. 

Figure  la  shows  a  schematic  of  a  current  sheet,  with  the  current  in  the  +x%  direction 
pointing  out  of  the  paper.  The  normal  component  Bz  =  Bn  =  const  of  the  magnetic  field 
is  in  the  direction  of  xz,  and  the  B\{xz)  component  in  the  X\  direction  reverses  direction 
at  xz  =  0.  We  set  Bz  —  0  in  the  present  study.  The  current  sheet  is  represented  by  the 
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region  between  the  dashed  lines  where  the  curvature  of  magnetic  field  “lines”  is  maximum. 
We  denote  the  half-width  of  this  region  by  S.  This  geometry  is  appropriate  for  a  distended 
dipole  field  such  as  the  magnetotail  and  the  heliospheric  current  sheet  or  for  a  rotational 
discontinuity.  We  will  use  subscripts  1,  2,  and  3  to  denote  vector  components  in  the  X\, 
X2,  and  £3  directions. 

In  the  solar-system  and  astrophysical  regimes  of  interest,  plasmas  are  nearly  collisionless, 
i.e.,  Vi  ft*.  Here,  Vi  and  ftj  are  the  characteristic  ion  collision  and  Larmor  frequencies, 
respectively.  More  specifically,  the  Larmor  radius  pn  =  t?th/ftn  of  thermal  ions  based  on 
Bn  and  the  time  ions  spend  in  the  current  sheet,  which  is  ~  ft”1,  are  much  shorter  than 
the  classical  mean  free  path  and  mean  free  time,  respectively,  where  vth  is  the  ion  thermal 
speed,  ft,»  =  eBn/miC,  mi  is  the  ion  mass,  and  c  is  the  speed  of  light.  In  these  regimes, 
the  structure  of  a  current  sheet,  whether  symmetric  or  asymmetric,  is  determined  by 
the  collisionless  motion  of  ions,  and  fluid  approximations  such  as  magnetohydrodynamics 
(MHD)  are  not  applicable. 

The  particle  motion  in  and  the  structure  of  a  current  sheet  can  be  approximated  as 
one  dimensional  (1-D)  provided  Li,Lv  »  pn,  where  Lx#  are  the  gradient  scale  lengths 
along  Xi  and  x2»  respectively  [Burkhart  and  Chen,  1993].  With  nonzero  Bn,  the  particle 
motion  has  only  two  constants  of  motion  in  involution  (the  Hamiltonian  and  one  canonical 
momentum)  and  is  nonintegrable  even  in  1-D  [Chen  and  Palmadesso,  1986].  The  charged- 
particle  orbits  consist  of  three  distinct  classes:  transient  (sometimes  called  Speiser-type), 
stochastic,  and  integrable  (regular)  orbits.  The  transient  and  stochastic  orbits  come  in 
from  and  escape  to  the  asymptotic  regions  outside  the  current  sheet,  defined  by  |a;3|  »  6, 
while  the  integrable  orbits  are  trapped  in  the  current  sheet.  In  Figure  la,  the  asymptotic 
source  regions  are  indicated  by  S+  (for  x3  »  +<5)  and  S~  (for  x3  -C  -S),  and  the  respective 
incoming  particle  distribution  functions  are  denoted  by  f+  and  /”.  The  time  spent  in  the 
current  sheet  is  ~  ft“1/2  for  transient  orbits  and  several  times  ft”1  for  stochastic  orbits. 
In  contrast,  the  motion  in  a  neutral  sheet  with  Bn  =  0  is  integrable  and  has  no  stochastic 
or  transient  orbits;  all  particles  in  such  a  neutral  sheet  remain  trapped.  This  is  true  with 
or  without  a  shear  component  B2  in  the  x2  direction. 

As  a  matter  of  notation,  we  will  use  subscript  a  to  denote  quantities  defined  in  the 
asymptotic  regions.  Accordingly,  Ba  =  Bi(|x3|  5)  is  the  asymptotic  value  of  the 

reversing  field.  For  current  sheets  of  interest,  Bn/ Ba  &  0.5.  Where  necessary,  we  will  use 
superscript  +  (— )  to  denote  quantities  defined  in  S+  (S~).  Where  the  distinction  is  not 
important,  we  will  use  S  to  denote  one  or  both  source  regions. 

The  theoretical  understanding  to  date  is  mostly  based  on  symmetric  current  sheets 
because  of  their  relative  simplicity  and  intrinsic  applicability  to  natural  systems  such  as 
the  magnetotail  that  are  topologically  symmetric.  In  an  early  theoretical  work,  Eastwood 
[1972]  studied  the  properties  of  1-D  symmetric  current  sheets  having  B2  =  0  and  Bn  ^  0, 
with  incoming  particle  distributions  (/+  =  /”)  consisting  of  cold  beams  in  which  all  par¬ 
ticles  have  identical  field-aligned  drift  speed  vd  at  the  asymptotic  sources.  The  treatment 
is  self-consistent  in  that  Ampere’s  law,  J  =  (c/47r)V  x  B,  is  satisfied.  This  work  was 
later  extended  to  include  thermal  spread  [Eastwood,  1974]  and  to  current  sheets  of  finite 
width  and  length  [Eastwood,  1975].  Francfort  and  Pellat  [1976]  analytically  studied  thin 
current  sheets  assuming  Bn/Ba  <  1  and  Vo/vth  »  1  and  using  an  adiabatic  invariant  of 
the  ion  motion  [Sonnerup,  1971].  In  these  studies,  only  the  transient  orbits  were  included. 
Lemaire  and  Burlaga  [1976]  constructed  steady-state  solutions  of  Vlasov-Maxwell  equar 
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tions  for  1-D  current  sheets  with  Bn  =  0  and  applied  them  to  tangential  discontinuities  in 
the  SW.  In  this  self-consistent  1-D  model,  the  current  sheet  is  allowed  to  be  asymmetric. 
With  Bn  =  0,  there  axe  three  constants  of  motion  in  involution,  and  all  trajectories  are 
integrable. 

More  recent  work  has  improved  on  the  above  models  by  including  the  full  particle 
dynamics  with  Bn  ^  0  using  simulation  techniques  [Burkhart  et  al.,  1992;  Pritchett  and 
Coroniti,  1992,  1995;  Holland  and  Chen ,  1993;  Harold  and  Chen ,  1996;  Hesse  et  al.,  1996] 
or  approximations  of  particle  orbits  assuming  Bn/Ba  ■<  1  [Kropotkin  et  al,  1997;  Sitnov 
et  al,  2000].  These  works  identified  the  quantity  vo/vt h  as  an  important  parameter, 
where  vth  is  the  thermal  speed  of  the  asymptotic  source  particle  distributions.  For  the 
magnetotail,  vo/vth.  1  is  appropriate  for  quiet- time  conditions,  while  vo/vth  1 
corresponds  to  field-aligned  source  distribution  functions  for  disturbed  situations  with 
strong  convection  electric  fields. 

All  of  the  above  equilibrium  models  produce  self-consistent  magnetic  field  profiles  sim¬ 
ilar  to  the  hyperbolic  tangent  Harris  field  [ Harris ,  1962],  but  the  plasma  properties  are 
noticeably  different  for  the  two  regimes.  In  the  vo/vth  1  regime  [Eastwood,  1972, 1974; 
Burkhart  et  al.,  1992;  Pritchett  and  Coroniti,  1992;  Kropotkin  et  al.,  1997],  the  particle 
density  is  highly  peaked  in  the  current  sheet,  as  in  the  Harris  model.  In  the  vo/vth  •C  1 
regime,  in  contrast,  the  particle  density  profile  is  nearly  constant  across  the  current  sheet 
[Holland  and  Chen,  1993;  Sitnov  et  al.,  2000].  The  latter  configuration  is  consistent  with 
the  observed  plasma  densities  in  the  quiet-time  magnetotail  current  sheet  [McComas  et 
al,  1986;  Sergeev  et  al.,  1993]. 

The  functional  form  of  the  incoming  source  distribution  is  also  an  important  parameter. 
Burkhart  et  al.  [1992]  used  a  drifting  Maxwellian  as  a  source  distribution  and  found  no 
equilibrium  solutions  for  their  simulations  in  the  vo/vth  &  1  regime.  In  contrast,  Holland 
and  Chen  [1993]  obtained  solutions  at  values  as  low  as  Vd/v^  =  0.05  using  a  ^-function 
distribution  [Christon  et  al,  1989],  which  is  prevalent  in  space  plasmas.  Physically,  the 
high-energy  tail  of  a  k  distribution  provides  excess  particles  executing  transient  orbits 
relative  to  a  Maxwellian  distribution,  producing  the  current  that  overcomes  plasma  dia¬ 
magnetism.  Sitnov  et  al.  [2000]  also  found  equilibrium  solutions  in  the  vo/vth  -C  1  regime 
assuming  Bn/Ba  -C  1.  In  these  studies,  the  pressure  tensor  was  found  to  be  anisotropic 
and  nondiagonal,  in  contrast  to  the  isotropic  and  diagonal  pressure  tensor  of  the  Harris 
model  while  Bi(x 3)  remained  similar  to  the  hyperbolic  tangent  Harris  field.  This  shows 
that  the  Harris-like  field  can  correspond  to  a  range  of  different  particle  distribution  func¬ 
tions  and  plasma  properties.  In  comparing  models  with  observations,  it  is  necessary  to 
examine  both  plasma  properties  and  the  magnetic  field. 

For  application  to  current  sheets  in  the  SW,  magnetopause,  solar  and  stellar  coronae, 
asymmetry  must  be  allowed.  For  the  SW,  two  common  types  of  discontinuities  have  been 
discussed  [e.g.,  Lepping  and  Behannon,  1986]:  rotational  discontinuities  (RDs),  character¬ 
ized  by  Bn  7^  0,  and  tangential  discontinuities  (TDs),  which  have  Bn  =  0  at  the  midplane 
of  discontinuity.  The  magnetic  topology  is  distinct  for  the  two  types  of  discontinuities, 
but  current  sheets  are  present  at  field  reversals  in  either  case.  Our  goal  here  is  not  to 
differentiate  between  the  different  types  of  magnetic  discontinuities,  but  rather  to  focus 
on  the  generic  properties  of  collisionless  current  sheets  at  magnetic  discontinuities. 

Lee  and  Kan  [1979]  used  a  method  similar  to  that  of  Lamaire  and  Burlaga  [1976]  and 
constructed  steady-state  solutions  of  the  Vlasov-Maxwell  equations  for  the  magnetopause 
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with  asymmetry  and  significant  jumps  in  physical  quantities.  Ding  et  al.  [1992]  investi¬ 
gated  the  dependence  of  the  collisionless  tearing  mode  on  the  asymmetry  of  the  current 
sheet.  Recently,  particle-in-cell  (PIC)  simulations  of  collisionless  magnetic  reconnection 
in  asymmetric  current  sheets  have  been  carried  out  [Swisdak  et  al.,  2003;  Pritchett,  2008]. 
Cassak  and  Shay  [2007]  provided  a  theory  and  simulation  of  asymmetric  reconnection  in 
the  MHD  regime.  Malova  et  al.  [2007]  proposed  a  model  of  asymmetric  thin  current 
sheets  using  approximate  ion  dynamics.  They  injected  an  ensemble  of  ions  from  one  side 
of  a  current  sheet  defined  by  B\,  B$  ^  0  and  B-j,  —  0  and  examined  the  spatial  structure 
of  Bi,  density  n,  and  current  J2  (all  expressed  in  our  coordinate  system  shown  in  Figure 
lb).  The  magnetic  field  profile  was  found  to  be  similar  to  the  hyperbolic  tangent  Harris 
field  for  varying  degrees  of  asymmetry.  The  self-consistency  of  the  results  and  the  ion 
flows  in  the  current  sheet  were  not  examined.  Mingalev  et  al.  [2009]  extended  this  work 
and  constructed  numerical  equilibria  satisfying  the  self-consistent  force  balance  required 
to  a  good  approximation.  In  both  models,  the  ion  density  profiles  have  pronounced  peaks 
at  the  field  reversal,  and  the  electron  contribution  to  the  current  is  neglected. 

In  the  present  paper,  we  develop  a  hybrid  kinetic  model  of  asymmetric  current  sheets 
including  the  electron  contributions,  treating  ions  as  kinetic  particles  and  electrons  as 
a  massless  fluid.  We  focus  on  the  detailed  particle  distributions  and  velocity  moments 
thereof.  The  model  current  sheet  is  self-consistent  and  is  characterized  by  an  anisotropic 
and  nondiagonal  pressure  tensor  as  well  as  asymmetric  flows.  The  calculated  current  sheet 
properties  will  be  directly  compared  with  high-resolution  solar  wind  data  from  the  WIND 
MFI  (0.044  sec  sampling  rate)  [Lepping  et  al,  1995]  and  WIND  3DP  (3-second  sampling 
rate)  [Lin  et  al,  1995]  measurements.  The  agreement  is  found  to  be  good.  The  apparent 
exception  is  the  ion  flow  speed  in  the  X2  direction. 

As  mentioned  earlier  in  connection  with  the  Harris  model,  the  magnetic  field  profile 
alone  is  not  sufficient  to  infer  current  sheet  properties.  The  SW  data  show  that  asymmetric 
and  sheared  particle  flows  across  current  sheets  are  correlated  with  observed  magnetic 
field  structures.  Thus,  any  model  that  attempts  to  explain  current  sheets  in  the  SW  must 
simultaneously  include  asymmetry  and  sheared  flows.  In  the  present  work,  we  extend  the 
self-consistent  test-particle  model  of  Holland  and  Chen  [1993]  by  including  asymmetry, 
strongly  sheared  flows  across  the  current  sheet,  and  an  approximate  treatment  of  the 
electron  contribution  to  the  current.  The  self-consistent  test-particle  ion  approach  is  the 
same  as  that  of  Mingalev  et  al.  [2009]  and  Holland  and  Chen  [1993]. 

The  remainder  of  this  paper  is  organized  as  follows.  In  section  2,  we  summarize  the 
basic  scales  governing  collisionless  current  sheets.  Section  3  discusses  the  current  sheet 
model,  with  the  model  results  given  in  section  4.  In  section  5,  we  compare  the  numerical 
results  with  WIND  spacecraft  data,  and  finally,  Section  6  gives  a  discussion  of  the  results 
along  with  our  conclusions. 

2.  Physics  of  Collisionless  Current  Sheets 

The  basic  physics  that  determines  the  structure  of  collisionless  current  sheets  is  the 
dynamics  of  the  charged  particles  in  the  current  sheet.  In  this  regard,  the  existence  of 
distinct  classes  of  orbits  corresponding  to  different  phase  space  regions  plays  an  essential 
role.  In  the  cold-beam  vd/vu,  >  1  regime  [Eastwood,  1972],  where  essentially  all  particles 
execute  identical  transient  orbits,  Francfori  and  Pellat  [1976]  obtained  the  scaling  for  the 
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where  pa  =  vth/$la  is  the  Larmor  radius  of  thermal  ions  in  the  asymptotic  field  component 
Ba,  and  fi0  =  eBJmiC.  This  expression  depends  on  an  approximate  adiabatic  invariant 
[Sonnerup,  1971]  for  transient  orbits  using  the  assumption  Bn/Ba  C  1  and  vo/vth  » 
1.  In  this  regime,  we  have  6  <  pa  (the  “zero-thickness”  limit).  The  scaling  was  verified 
for  vD/vth  ^  1  by  a  numerical  simulation  [Burkhart  et  al.,  1992].  In  contrast,  in  the 
vo/vth  <  1  regime,  a  heuristic  derivation  [Chen  et  al,  1990]  yielded  the  scaling 

S  ~  pa,  (2) 

in  agreement  with  the  results  of  a  more  detailed  analysis  including  the  nonintegrable 
motion  [Sitnov  et  al.,  2000]  and  the  test-particle  simulation  of  Holland  and  Chen  [1993]. 
These  scaling  properties  are  universal  in  1-D  in  that  they  only  depend  on  the  collisionless 
ion  motion  in  a  current  sheet.  Note  that  pa  can  be  related  to  the  ion  skin  depth  A,  in  the 
asymptotic  region  by 

(3) 

where  A*  =  c/ojpi,  ujp,  =  (47rn,e2/m<)1^2  is  the  ion  plasma  frequency,  and  0a  = 
2nT/(B%/8n),  assuming  Te  =  T,  =  T  for  electron  and  ion  temperatures.  Here,  the 
Boltzmann  constant  has  been  absorbed  into  T,  and  we  will  express  temperature  in  units 
of  energy.  For  0a  =  1,  we  have  pa  =  \i/y/2. 

The  different  vo/vth  regimes  are  determined  by  the  asymptotic  source  distribution  func¬ 
tions,  which  control  the  relative  population  of  transient  and  stochastic  particles.  Although 
the  particle  motion  cannot  be  directly  observed,  certain  observed  features  such  as  peaks 
and  valleys  in  ion  velocity  distributions  obtained  in  the  quiet-time  magnetotail  [Huang 
et  al.,  1989;  Holland  et  al.,  1999]  can  be  accurately  explained  as  a  direct  manifestation 
of  the  characteristic  phase  space  resonance  exhibited  by  transient  and  stochastic  orbits 
[Chen  et  al.,  1990;  Burkhart  and  Chen,  1991;  Wang,  1994;  Savenkov  and  Zelenyi,  1996; 
Holland  et  al.,  1999].  This  suggests  that  the  assumption  of  collisionless  particle  motion 
is  valid  in  the  magnetotail.  The  physics  scales  with  the  ion  Larmor  radius,  implying  that 
the  structure  of  collisionless  current  sheets  can  be  scaled  from  the  magnetosphere  to  the 
solar  wind  as  well  as  to  laboratory  plasmas,  where  it  can  be  tested  against  directly  mea¬ 
surable  plasma  and  magnetic  field  data,  and  ultimately  to  the  solar/stellar  coronae.  It  is 
the  purpose  of  the  present  work  to  (1)  understand  the  structure  and  (2)  test  the  physics 
of  thin  collisionless  current  sheets  in  the  solar  wind. 

Note  that  the  usual  plasma  0  =  2nT/(B2/8n)  is  based  on  the  total  magnetic  field  and 
is  of  order  unity  or  greater  in  a  current  sheet  even  if  the  medium  has  typical  0  1, 

provided  there  is  no  strong  “applied”  magnetic  field  due  to  current  distributions  outside 
the  current  sheet. 

3.  Asymmetric  Current  Sheet  Model 

The  model  obtains  numerical  solutions  to  the  steady-state  Vlasov-Maxwell  equations 
satisfying  J  =  (c/47t)VxB.  In  this  approach,  the  orbits  of  ions  are  obtained  by  integrating 
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the  full  equations  of  motion  in  electric  and  magnetic  fields  and  summing  the  contribution 
from  each  ion  according  to  specified  input  distributions,  /+  and  /“.  The  model  is  one- 
dimensional,  having  spatial  variation  only  in  the  direction  perpendicular  to  the  plane  of 
the  current  sheet. 

3.1.  Coordinate  System 

Figure  lb  shows  the  coordinate  systems  used  in  the  analysis.  The  box  represents  the 
current  sheet  and  corresponds  to  the  region  of  thickness  26  between  the  dashed  lines  in 
Figure  la.  The  x-y-z  coordinates  are  centered  at  the  observer  (e.g.,  the  WIND  spacecraft) 
such  that  the  unit  vectors  x,  y,  and  z  are  aligned  with  those  of  the  usual  Geocentric  Sun- 
Earth  (GSE)  coordinates.  In  this  frame,  +x  points  from  the  Earth  to  the  Sun,  +y  points 
out  of  the  page,  and  +z  points  northward.  The  xi-x2-X3  coordinate  system  is  defined  to 
be  co-moving  with  the  current  sheet  with  the  X3  axis  normal  to  the  plane  of  the  current 
sheet  and  the  x\  axis  parallel  to  the  reversing  component  of  the  magnetic  field.  The  dots 
on  the  Xi,  x2,  and  x3  axes  axe  intended  to  help  visualize  where  the  three  axes  pierce  the 
surfaces  of  the  box  facing  out  of  the  page  at  various  angles. 

Figure  lb  depicts  the  current  sheet  (thickness  26)  as  moving  in  the  —x  direction  (to  the 
right)  in  the  x-y-z  frame  with  a  velocity  Vsw-  If  the  structure  intercepts  the  observer  on 
the  x  axis,  the  asterisk  (*)  to  the  right  represents  the  entry  point  in  the  lower  surface  of 
the  box  and  the  left  asterisk  the  exit  point  in  the  upper  surface  of  the  box.  The  distance 
between  the  asterisks  is  always  greater  than  or  equal  to  the  actual  thickness  26  and  can 
be  much  greater  depending  on  angle  of  inclination  of  the  current  sheet  to  the  x-y  plane. 

The  two  coordinate  systems  are  related  by 

r  =  72(x  +  xr)  (4) 

where  x  =  (x, y, z),  r  =  (xi,  x2, X3),  xT  accounts  for  the  translational  motion  of  the 
current  sheet  with  respect  to  the  observer,  and  72  is  a  matrix  specifying  rotations  in  the 
observer  {x-y-z)  coordinate  system.  The  vector  xT  is  necessary  because  the  observed  data 
are  obtained  as  a  time  series  as  the  structure  advects  past  the  observer  at  velocity  V$vv- 
Thus,  we  convert  the  time  series  to  a  function  of  x  by  using 

xr  =  —V  swt, 

where  t  is  time  and  Vsvv  is  a  constant  velocity  vector  chosen  to  be  the  SW  velocity  at 
the  center  of  the  current  sheet.  The  transformation  matrix  is  standard  and  is  given  by 
three  successive  rotations  [Goldstein,  1980] 

72.  =  'R-x"  {6x" )  72z/  (#z< )  IZy  (#y )  ■  (5) 

Here,  7 Zy  is  the  rotation  through  angle  9V  about  the  GSE  y  axis  centered  at  the  observer. 
The  resulting  coordinates  are  denoted  by  x7,  y',  and  z! .  The  operator  TZz>  is  a  rotation 
through  0Z>  about  the  z'  axis.  The  resulting  coordinates  are  x",  y",  and  z".  Similarly, 
72j»  performs  a  rotation  through  angle  6X»  about  the  x77  axis.  This  yields  the  desired 
coordinate  system,  which  is  denoted  by  xi,  x2,  and  X3.  We  carry  out  the  simulations  in 
the  xi-x2-X3  coordinates.  The  observational  data  are  transformed  from  the  GSE  into  the 
these  coordinates  for  comparison  with  the  model  results. 
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3.2.  Self-Consistent  Equilibrium 

The  model  structure  is  an  equilibrium  current  sheet  that  is  stationary  in  the  Xi-X2-x3 
coordinate  system.  For  simplicity,  we  assume  that  the  fields  are  translationally  invariant 
in  the  X\  and  X2  coordinates,  having  variation  only  in  the  normal  (x3)  direction.  This  1-D 
approximation  is  valid  provided  Li>2  pn,  where  L\^  are  the  gradient  scale  lengths  in 
the  Xi  and  x2  directions,  respectively. 

The  equilibrium  structure  of  the  current  sheet  is  determined  by  the  motion  of  ions.  The 
motion  and  therefore  the  equilibrium  Vlasov  distribution  functions  axe  functions  of  the 
constants  of  motion.  Although  the  magnetic  field  is  translationally  invariant  in  xi  and  x2, 
the  presence  of  nonzero  B3  =  B„  means  that  the  Hamiltonian  necessarily  depends  on  x\ 
or  x2  in  addition  to  x3,  rendering  the  particle  motion  nonintegrable  and  chaotic  (section 
1).  We  choose  a  gauge  in  which  the  canonical  momentum  p2  =  niiV2  +  qiA2{x\,x2)/c 
conjugate  to  X2  is  conserved,  where  A  is  the  vector  potential. 

To  start  the  simulation  an  initial  magnetic  field  B  and  an  electric  field  E  are  speci¬ 
fied.  For  simplicity,  we  assume  that  the  magnetic  field  is  unsheared,  i.e.,  B2  =  0.  In  the 
Vohth  >  1  regime,  hybrid  simulations  [e.g.,  Pritchett  and  Coroniti,  1992]  indicate  that 
a  nonzero  B2  component  in  the  form  of  radiating  waves  is  generated.  In  the  Vp/vth  < 
1  regime,  in  contrast,  any  B2  component  that  may  arise  is  small  and  grows  slowly  (hun¬ 
dreds  of  fi"1)  [Cargill  et  al.,  1994],  so  that  the  assumption  of  unsheared  magnetic  field  is 
consistent  with  uoA'th  1.  Thus,  we  take  the  magnetic  field  to  have  the  form 

B(x3)  =  Bi(x3)xx  +  B„x3,  (6) 


where  Bn  is  constant  and  x*  denotes  the  unit  vector  in  the  x*  direction.  With  B2  =  const 
(=  0)  and  d/dxi  =  d/dx2  =  0,  the  total  current  density  J  is  entirely  in  the  X2  direction, 


J=fVxB  = 
47T 


c  dBi^ 
47t  dxz  Xa’ 


(7) 


with  J2{xz)  =  (c/4n)dBi/dxz  =  Jziiz^+Jieixz)-  In  the  hybrid  approach,  Jn  is  calculated 
using  individual  ion  orbits  in  order  to  include  the  important  Larmor  radius-scale  kinetic 
effects,  while  J2e  is  calculated  using  the  electron  fluid  equations,  taking  electron  inertia 
to  be  negligible  on  the  ion  time  scales  ~  fi"1. 

In  a  1-D  equilibrium,  Vx  E  =  0  implies  that  E\  and  E2  are  both  constant.  The 
component  Ez  is  determined  by  Poisson’s  equation 

V-E  =  =  4tt  eSn,  (8) 

OXz 

where  Sn  =  n*  —  ne.  Because  Ez  varies  over  the  current  sheet  with  thickness  of  the 
order  of  p0,  it  can  be  shown  (Appendix  A)  that  Sn/n  ~  (Vk/c)2  1,  where  Vk  is  the 

characteristic  Alfv6n  speed.  Thus,  charge  neutrality  is  valid  to  lowest  order,  and  with  this 
understanding  we  use 


m  =  ne  =  n.  (9) 

Note  that  charge  neutrality  is  not  an  explicit  assumption  in  this  model.  Rather  it  depends 
on  Va/c  1  and  the  current  sheet  thickness  being  pa  as  determined  by  the  ion  motion. 
In  Lemaire  and  Burlaga's  model  with  Bn  =  0,  the  scalar  potential  and  charge  separar 
tion  are  included,  finding  charge  separation  to  be  extremely  small  ( Sn/n  ~  10-8)  in  the 
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Vjj/vth  1  regime.  Lee  and  Kan  [1979]  assumed  quasi-neutrality,  neglecting  ( V/c )2  -C  1 
contributions.  In  our  model,  the  scalar  potential  is  neglected. 

We  use  an  iterative  method  to  compute  the  self-consistent  magnetic  field,  Bi,  and 
the  corresponding  plasma  properties  of  the  current  sheet.  In  this  algorithm,  E2  and 
Bn  are  held  fixed,  but  Bifa)  is  iterated  to  self-consistency  satisfying  Ampere’s  law, 
J  =  (c/4tt)  V  x  B.  A  uniform  E2  field  can  be  transformed  away  by  a  so-called  de  Hoffman- 
Teller  transformation  in  the  x%  direction  so  that  the  calculation  can  be  done  by  setting 
E2  —  0.  The  de  Hoffman-Teller  transformation  velocity  is  Vjt  =  cE2/Bn,  which  then 
must  be  added  to  the  calculated  V2  if  E2  ±  0.  We  will  keep  E2  in  the  formulation  for 
possible  generalization  to  the  nonuniform  E2  ^  0  case,  but  the  actual  simulation  will 
be  done  using  E^  =  0  everywhere.  The  component  E\,  which  is  not  affected  by  the  de 
Hoffman-Teller  transformation,  is  set  to  zero  because  of  the  high  electrical  conductivity 
of  a  collisionless  plasma.  In  the  presence  of  asymmetry,  there  is  generally  a  net  pressure 
variation  across  the  current  sheet,  so  that  there  can  be  a  nonzero  E3  component  in  the 
current  sheet  if  the  pressure  tensor  is  symmetric.  If  the  pressure  tensor  is  allowed  to  be 
nondiagonal,  however,  it  is  possible  to  construct  self-consistent  equilibria  with  E3  =  0. 
In  the  present  paper,  we  limit  our  consideration  to  a  subclass  of  equilibria  having  E3  — 
0.  A  number  a  general  properties  and  implications  pertaining  to  this  simplification  are 
discussed  in  section  3.5  and  Appendix  A.  Note,  however,  that  the  treatment  in  Appendix 
A  is  heuristic.  In  the  following  sections,  we  will  discuss  in  detail  our  treatment  for  the 
ions,  electrons,  and  the  computation  of  the  self-consistent  fields. 


3.3.  Contribution  of  Ions  to  Equilibrium 

The  ions  are  treated  as  particles.  For  each  iteration  of  the  equilibrium  computation, 
a  specified  ensemble  of  ions  is  injected  from  the  boundaries  of  the  simulation,  and  the 
particle  trajectories  are  calculated  using  the  equations  of  motion 


v  = 


dX 

dt 


dv 

~dt 


§■  [e(X3)  +  -  x  B(X3) 

77lj  L  c 


(10) 


where  *  indicates  the  ion  species,  is  the  ion  charge,  mi  is  the  ion  mass,  and  c  is  the  speed 
of  light.  Here  and  henceforth,  we  use  X(f)  and  v(t)  to  denote  the  position  and  velocity  of 
ions  along  their  trajectories,  where  t  is  time  (an  orbit  parameter).  The  B  and  E  fields  are 
held  fixed  in  time  during  each  complete  iteration  of  the  equilibrium  computation,  which 
starts  when  the  particles  are  injected  into  the  system,  and  ends  when  all  the  particles 
have  exited  from  the  system. 

For  each  iteration,  the  contribution  of  each  particle  to  the  first  three  velocity  moments, 
density  ni(x3),  current  density  J<(x3),  and  the  energy  density  tensor  Qi(x$),  is  computed 
using  the  particle  weighting  function  iot(vtn)  for  each  particle.  The  weighting  function 
includes  both  the  values  of  the  asymptotic  source  density  (nf)  and  a  velocity  space  factor 
which  is  a  function  of  the  particle’s  initial  velocity,  v<n,  when  it  is  injected  into  the  system. 

The  particle  orbits  and  the  velocity  moments  are  computed  on  a  finite  grid  rather  than 
being  continuous  functions,  and  linear  interpolation  is  used  to  distribute  particle  quantities 
onto  the  grid.  For  this  purpose,  we  use  the  coordinates  Xi,x2,  and  x3  to  represent  a  finite, 
discrete  set  of  grid  coordinates  with  uniform  grid  spacing  of  Ax.  (Note  that  the  equations 
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of  motion  are  not  discretized.)  Where  multiple  species  of  ions  axe  included,  each  ion 
species  is  represented  by  M<  particles  in  the  equilibrium  simulation,  and  the  jth  particle 
is  represented  by  N&  discrete  orbital  points.  For  the  jth  particle,  these  orbital  points 
are  sampled  at  discrete  times  tm  and  denoted  by  X3(tm),  where  m  =  1,2,  The 

sampling  time  interval,  At  =  tm+i  -  fm,  is  equally  spaced.  The  total  contribution  to  the 
equilibrium  velocity  moments  at  a  grid  point  £3  is  obtained  by  summing  the  contributions 
from  all  the  particles  that  pass  within  a  grid  spacing  Ax  of  £3.  The  contribution  of  the 
jth  particle  of  species  i  to  the  equilibrium  velocity  moments  is  given  by 

1  nU) 

niU)(x3)  =  Y  «'i(i)(Vin)dv<nC(£3,  X3(tm)), 

•*v  m=l 

1  N<J) 

Ji0)(£3)  =  Jjyj  Y  9<vwp')(vin)dv<nC(£3,  X3{tm)),  (11) 

m=l 

1  N(J)  ,.v 

Q<0)(£3)  =  JjJjj  Y  "kVVWi  \vin)dvinC{x3,X3(tm)), 

m=l 

where  dvin  is  the  differential  volume  element  in  velocity  space  and  Qj  is  a  tensor.  The 
scalar  function  C(£3,  X3 (tm))  is  used  to  interpolate  between  the  particle  positions  and  the 
grid  cell  coordinates,  and  is  given  by 


C(x3,  X3(tm))  —  0,  (£3  X3\  >  Ax, 

C(x 3,  Xafc,))  =  |x3  -  X3|  <  fix. 

The  equilibrium  profiles  for  the  first  three  moments  for  each  ion  species  are  then  calculated 
by  summing  over  all  the  particle  trajectories  for  a  particular  species.  Thus,  the  number 
density  for  the  ith  species  is  given  by: 


1  M> 

(xa)  = 

mi  j= 1 

where  Mi  is  the  number  of  particles  representing  species  i.  The  total  ion  current  density 
Ji(x3)  and  energy  density  tensor  Qi(x3)  are  obtained  by  analogous  summations  over  j. 

From  these  velocity  moments,  other  quantities  such  as  the  equilibrium  flow  velocity,  Vj 
and  the  equilibrium  pressure  tensor  Pi  can  be  determined  as  follows: 

V  =  J-j 

*  qiXii  *’ 

Pi  =  Qi  -  miTiiViVi.  (12) 

The  ion  temperature  is  then  computed  from  the  trace  of  the  ion  pressure  tensor 

tr{  Pi} 


Ti  = 


3n i 


(13) 


In  the  equilibrium  current  sheet  the  force  balance  for  each  ion  species  is  given  by 

V  •  (  Pi  +  miTiiViVi)  -  ntfi  (e  +  £ V<  x  b)  =  0,  (14) 
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where  Pi  is  the  ion  pressure  tensor,  V*  is  the  average  velocity,  and  rii  is  the  ion  number 
density.  The  Lorentz  force  (E+V*  x  B/c)  is  balanced  by  kinetic  forces  V  •  (  Pi+mjniVjV<) 
for  each  species  independently  since  the  effect  of  collisions  between  ions  of  different  species 
are  small  on  the  relevant  time  scale  and  are  neglected. 

For  a  system  with  multiple  ion  species,  all  ion  contributions  are  summed  over  i.  Hence¬ 
forth,  however,  we  will  include  only  one  positive  ion  species,  i.e.,  protons,  and  designate 
all  proton  quantities  by  subscript  i. 

3.4.  Contribution  of  Electrons  to  Equilibrium 

We  expect  the  ion  dynamics  to  determine  the  basic  length  scale  in  the  problem,  i.e., 
the  thickness  of  the  current  sheet,  and  for  this  reason  (and  because  it  would  require 
much  greater  computing  resources)  we  do  not  include  the  individual  electron  orbits  in 
the  model.  The  model  does,  however,  include  the  electron  contribution  to  the  net  (total) 
current,  treating  the  electrons  as  a  massless  fluid.  The  electron  fluid  momentum  equation 
is  given  by: 

neTnelt  =  0  =  ~eHe  (E  +  \Ve  x  B)  -  v  •  Pe  (15) 

where  Vc  is  the  electron  flow  velocity,  Pe  is  the  electron  pressure  tensor,  and  the  electron 
mass  me  is  set  to  zero.  Our  treatment  of  the  electrons  is  similar  to  the  hybrid  simulation 
of  Leroy  et  al.  [1982],  except  that  their  simulation  is  time  dependent  while  our  calculation 
yields  time-independent  solutions. 

From  Ampere’s  law  (7),  the  1-D  system  satisfies  J3(x3)  =  0  everywhere.  With  condition 
(9),  this  gives  «/3  =  en(V3i  —  14.)  =  0  and,  therefore,  V&  =  V&.  =  V3.  Equation  (7)  also 
shows  that  J\  (x3)  =  0  for  B2  =  const  regardless  of  the  constant  value,  with  the  result 
Vu  =  Vie  =  Vi-  The  total  current  density  including  the  electron  contribution  is 

J2(x3)  =  en(V2i  -  V2e).  (16) 

The  electron  contribution  to  the  net  current  is 

J<le  —  CTiV^g, 

where  is  determined  by  the  electron  momentum  equation  (15)  and  n  =  ne. 

The  X3  component  of  (15)  is 

-en,(i*-iv*Bi)-^  =  0,  (17) 

and  the  x3  component  of  the  steady-state  ion  momentum  equation  is  obtained  from  equa¬ 
tions  (12)  and  (14): 

e^-iv*Bl)-^S  =  0.  (18) 

Eliminating  E3  from  equations  (17)  and  (18)  and  using  equation  (9),  we  obtain 

C  /  dQtti  d^33e\ 

enBi  \  8x3  8x3  ) 


2e  — 


(19) 
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For  the  electron  fluid,  we  use  the  polytropic  equation  of  state  [Krall  and  Trivelpiece, 
1986] 


where  p  =  nT  is  the  scalar  fluid  pressure,  p  =  mn  is  the  mass  density,  and  7  is  the 
polytropic  index.  Assuming  that  our  electron  fluid  obeys  a  similar  equation  of  state,  we 
obtain 


(20) 


where  7e  refers  to  the  electron  fluid.  Here,  we  have  assumed  diagonally  dominant  and 
nearly  isotropic  electron  pressure  tensor  so  that  =  pe.  This  is  reasonable  on  the 
spatial  and  temporal  scales  of  ion  dynamics  because  of  the  much  faster  electron  time 
scales.  In  addition,  since  we  assume  steady-state  (d/dt  =  0),  we  have  d/dt  =  V^d/dx^, 
and  equation  (20)  becomes 


(21) 


Using  equations  (16),  (19),  and  (21),  we  obtain  the  total  electric  current 


where  we  have  used 


(22) 


Here,  we  recognize  Joi  as  the  diamagnetic  current  determined  by  the  ion  motion,  and  the 
second  term  in  the  square  brackets  is  the  electron  contribution  to  the  total  current.  In 
obtaining  equation  (22),  we  have  made  the  approximation  Pe/Pz&i  ~  Te/Tj.  Thus,  the 
total  current  J2  can  be  determined  from  the  ion  quantities  if  %  and  Te/Ti  are  given. 

In  the  solar  wind,  the  density  scale  length  across  the  current  sheet  tends  to  be  compara¬ 
ble  to  the  pressure  scale  length,  and  the  electron  and  ion  temperatures  axe  also  comparable. 
The  (Te/Ti)  factor  in  equation  (22)  shows  that  the  electron  contribution  to  the  current 
can  be  significant.  In  the  quiet-time  magnetotail  current  sheet,  on  the  other  hand,  we 
have  Te/T<  <C  1,  and  the  electron  contribution  is  less  important. 


3.5.  The  Self-consistent  Equilibrium  Fields 

For  each  iteration,  the  total  equilibrium  current  is  computed  using  the  ion  current 
JDi  and  equation  (22).  In  order  to  obtain  a  smoother  solution  and  improve  the  rate  of 
convergence  in  the  numerical  algorithm,  we  average  the  gradients  in  equation  (22)  over 
the  region  containing  the  current  sheet,  yielding 


J  =  Jm 


hmmmi) 


*2, 


(23) 


where  (•  •  •)„  indicates  the  average  over  the  current  sheet.  We  assume  an  adiabatic  electron 
fluid  and  take  7e  =  5/3.  In  the  actual  simulation,  we  use  Joi  =  enV ru  to  improve  numerical 
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convergence.  This  substitution  results  in  a  local  error  of  Joi  —  Ja  =  —cen(E3/Bi).  The 
E3  component  can  be  obtained  from  the  ion  momentum  (Appendix  A),  yielding 


E3  = 


1 

enBn 


dQm  D  dQ33i 

D\  +  — — tf3 


dx3 


dx3 


(24) 


This  shows  that  E3,  which  is  parallel  to  Bn,  depends  on  pressure  gradients.  In  Harris-like 
equilibria,  the  off-diagonal  tensor  elements  are  zero,  and  an  electric  field  E3  is  necessary 
to  balance  the  forces  in  the  x3  direction.  In  our  equilibrium  system,  where  the  pressure 
tensor  is  allowed  to  be  nondiagonal,  force  balance  can  be  achieved  with  the  gradients  of 
off-diagonal  tensor  elements.  In  our  simulation,  we  specialize  to  a  class  of  self-consistent 
equilibria  with  E3  =  0.  In  such  a  system,  ( dQ3u/dx3)Bi  balances  (dQ^/dx^B^  In 
this  regard,  we  note  that  numerically  computed  E3,  with  derivatives  of  Q3ki,  tend  to 
be  numerically  noisy  so  that  setting  £3=0  allows  rapid  convergence.  We  also  note 
that  averaged  over  the  current  sheet,  the  two  extrema  in  E3  nearly  cancel,  and  the  net 
contribution  of  E3  to  the  total  current  calculation  mainly  depends  on  the  degree  of  overall 
asymmetry  in  the  pressure  across  the  current  sheet  (Appendix  A).  This  is  expected  to  be 
small  in  our  case. 

Once  computed,  the  total  J  is  used  to  calculate  the  magnetic  field  B  from  Ampere’s 
law  V  x  B  =  (47t/c)J.  The  newly  calculated  B  is  then  used  as  the  new  prescribed 
magnetic  field  for  the  computation  of  particle  orbits.  The  procedure  is  iterated  until 
the  difference  between  magnetic  fields  of  successive  iterations  is  negligible.  In  practice, 
we  use  a  relaxation  technique  in  which  we  mix  the  old  and  newly  calculated  magnetic 
fields  in  order  to  generate  a  smoother  result  and  minimize  overshooting  from  iteration  to 
iteration.  Thus,  the  magnetic  field  B  for  the  equilibrium  current  sheet  is  self-consistently 
determined  for  an  imposed  E2  electric  field  and  particle  velocity  distributions  specified  in 
the  boundary  regions. 

Heretofore,  we  have  formally  carried  E  to  show  the  general  structure  of  the  equations. 
In  the  actual  calculations,  however,  we  assume  that  E  =  E2x 2  with  E2  —  const  and  apply 
a  de  Hoffman-Teller  transformation.  That  is,  the  calculations  will  be  performed  in  the 
frame  in  which  E2  =  0  everywhere.  In  the  original  frame,  an  E  x  B  translation  velocity, 
Vjr  =  cE2/Bn  is  added  to  the  velocity  of  each  particle. 


4.  Numerical  Results 

For  the  initial  guess  of  Bi(x3),  we  use  equation  (6)  with  the  Harris  profile, 

Bi  (x3)  =  Bao  tanh(x3/ tf0)  (25) 

where  B^  is  the  initial  value  of  the  asymptotic  B\  field,  <J0  is  the  initial  thickness  of  the 
magnetic  discontinuity  corresponding  to  the  initial  half-width  of  the  current  sheet. 

The  Bi(x3)  component  of  the  magnetic  field  evolves  with  each  iteration  of  the  calculation 
until  convergence  of  the  magnetic  field  is  reached.  The  scale  length  of  the  converged 
magnetic  field  discontinuity,  B\/(dBi/dx3),  can  be  quite  different  from  60 •  We  have 
found  the  converged  magnetic  field  B\ (x3)  to  be  insensitive  to  both  the  initial  Boq  and  So- 
The  simulation  encompasses  three  regions.  The  current  sheet  is  centered  at  x3  =  0  and 
is  located  in  the  central  region  given  by  —Lc  <x3<  Lc,  where  Lc  »  <50.  We  typically  use 
Lc  which  is  10-20  8q.  On  either  side  of  the  central  region  is  an  asymptotic  source  region, 
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defined  by  Lc<x 3<  La  (S+)  and  —La  <  x3  <  —Lc  (S~).  Thus,  the  actual  simulation  is 
bounded  by  the  outer  surfaces  of  the  asymptotic  regions,  x3  —  ±La.  For  the  asymptotic 
regions,  we  choose  (La  —  Le)  >  3pa. 

In  the  simulation,  we  prescribe  two  ion  distribution  functions,  /+  and  f~,  and  inject 
particles  accordingly  at  the  center  of  the  S+  and  S~  regions.  They  are  then  allowed  to 
move  into  the  centred  region,  —Lc  <  x3  <  Lc.  The  asymptotic  regions  S  are  provided  to 
allow  the  particles  to  follow  their  correct  dynamics  beyond  the  limits  (£3  =  ±LC)  of  the 
central  region  until  they  eat  the  simulation  through  one  of  the  boundaries  at  x3  =  ±La. 
Because  we  inject  particles  from  the  asymptotic  regions,  there  are  no  integrable  orbits 
(trapped  particles)  except  for  a  small  number  that  may  be  scattered  into  the  integrable 
regions  [Holland  and  Chen,  1993]  in  the  numerical  computation.  A  quantitative  analysis 
of  the  effect  of  trapped  particles,  which  is  not  expected  to  be  significant,  is  beyond  the 
scope  of  this  paper. 

For  the  incoming  source  distribution  functions  in  the  S+  (S~)  asymptotic  region,  the 
average  flow  speed  is  specified  to  be  (f/f )  in  the  xi  direction.  We  use  a  Ac-like  initial 
velocity  space  distribution  [Christon  et  al.,  1989],  consistent  with  observations  of  SW 
plasmas  [e.g.,  Gloeckler  et  al,  1992;  Chotoo  et  al.,  1998].  A  ^-function  is  a  Maxwellian  at 
low  energies,  having  a  power  law  tail  at  high  energies.  Defining  v2  =  (vi  —  U\)2  +  v\  +  v\, 
we  use  a  «-like  distribution  function  defined  by 


f(v)  =  no  exp 


v  <  2vth, 


f(v)  =  no 


v  >  2 vth, 


(26) 


where  no  is  the  density  normalization  factor  and  vth  is  the  thermal  speed  of  the  Maxwellian. 
In  our  study,  we  will  use  k  =  5.  It  is  understood  that  U%,  no,  and  vth  are  separately  defined 
for  S+  and  S~.  Note  that  Ui  is  roughly  analogous  to  the  parameter  Vp  of  the  symmetric 
case,  which  is  essentially  the  specified  field-aligned  drift  velocity  in  the  asymptotic  region. 
For  the  asymmetric  case,  however,  we  have  U+  /  C/f ,  and  there  is  no  simple  parameter 
that  precisely  corresponds  to  vp.  In  some  cases,  such  as  the  example  to  be  considered 
here,  the  average  of  U+  and  t/f  can  roughly  play  the  role  of  vp. 

The  particle  distribution  that  can  be  measured  by  an  observer  is  the  total  distribu¬ 
tion  consisting  of  the  input  and  output  distributions,  the  latter  being  determined  by  the 
numerical  integration  of  orbits.  We  define  Vj+  (V^)  as  the  average  ion  velocity  in  the 
X\  direction  of  the  total  particle  distribution  in  asymptotic  region  S+  (S~).  The  total 
particle  distribution  functions  in  S±  also  yield  other  velocity  moments  such  as  the  ion 
density  nf  and  temperature  T±. 

Let  us  now  consider  an  illustrative  current  sheet  computation  appropriate  for  a  SW 
plasma  near  1  AU.  Figure  2  describes  a  converged  equilibrium  current  sheet  obtained  by 
our  iterative  method.  Anticipating  comparison  of  numerical  results  with  SW  data,  this 
solution  has  been  obtained  to  yield  the  best  overall  fit  to  a  SW  current  sheet  to  be  discussed 
in  the  next  section.  From  the  top,  the  figure  shows  the  reversing  component  £1(0:3),  the 
total  (ions  plus  electrons)  current  density  J3{x3),  ion  density  n<(x 3),  and  ion  temperature 
Ti(x 3)  in  the  central  simulation  region,  —  Lc  <  x3  <  Lc,  where  we  have  chosen  Lc  =  4000 
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km.  The  Bi  profile  can  be  approximately  by  a  hyperbolic  tangent  function  and  resembles 
the  Harris  field.  The  half-width  thickness  of  converged  the  current  sheet  is  S  ~  350  km, 
which  is  approximately  three  times  the  ion  thermal  Larmor  radius  =  vth/^a  —  120  km 
at  X3  =  La.  This  confirms  the  scaling  relation  (2)  and  the  importance  of  using  a  kinetic 
approach  for  this  problem.  The  current  density  J2(x 3),  plotted  in  the  second  panel,  is 
highly  localized  within  the  central  region  and  shows  a  characteristic  peak  at  the  midplane, 
with  J2  ~  0  elsewhere.  The  third  panel  shows  ion  density,  Ui(x 3),  which  exhibits  an 
enhancement  of  about  25%  in  the  current  sheet  in  comparison  with  the  asymptotic  value, 
n+  ~  3  cm-3,  of  the  converged  solution.  The  ion  temperature,  given  by  equation  (13),  is 
plotted  in  the  bottom  panel.  Of  particular  note  is  the  anti-correlated  asymmetry  that  is 
present  in  the  temperature  and  density  profiles. 

Note  that  the  ion  skin  depth  A *  for  this  example  is  comparable  to  p0.  Using  the  asymp¬ 
totic  density  of  the  converged  solution,  n+  ^  3  cm-3,  we  obtain  A <  ~  130  km,  showing 
pa  ~  A i.  The  asymptotic  values  ~  3.7  nT  and  ~  10  eV  yield  /?+  ~  2,  consistent 
with  equation  (3).  The  usual  plasma  /?  is  defined  with  the  total  magnetic  field.  For  the 
parameters  given  above,  /3+  =  2nT/(B2/8ir)  ~  1.6  in  region  S+  for  the  converged  current 
sheet.  At  the  center  of  the  current  sheet  where  B\  =  0,  we  have  local  (3  «  40. 

For  this  example,  the  quantity  (•  •  •)<*  in  equation  (23)  is  computed  to  be  ~  0.6.  Using 
7e  =  5/3  and  Te  =  Ti,  we  find  that  J2e  —  Ja  so  that  the  electron  current  accounts  for 
approximately  half  of  the  total  current  J2.  In  the  magnetopause  model  of  Lee  and  Kan 
[1979]  where  Ti/Te  ~  10,  the  electron  current  was  found  to  be  about  one  quarter  of  the 
total  current. 

The  asymmetry  in  n,  (x 3)  and  Tj(x 3)  arises  primarily  from  the  asymmetric  input  source 
particle  densities.  Force  balance  may  also  require  asymmetry  in  ni(x 3)  and  Ti(x3)  if  there  is 
asymmetry  in  the  flow  velocity.  That  is,  if  the  directed  momentum  flow  (njmiViVj  term) 
is  not  balanced  across  the  current  sheet,  an  asymmetry  in  the  temperature  is  required  to 
achieve  total  force  balance. 

Figure  3  shows  the  components  of  the  average  ion  flow  velocity  Vi,  obtained  as  the 
first-order  velocity  moment  of  the  total  ion  distribution  function  in  the  simulation  current 
sheet.  The  Vh(x3)  profile  is  asymmetric  and  nonzero.  In  contrast,  in  the  symmetric  case, 
everywhere  [e.g,  Holland  and  Chen,  1993;  Harold  and  Chen,  1996].  The  asymmetry 
in  the  average  velocity  is  mainly  due  to  Ui  ^  t/f.  The  second  panel  shows  that  V2 i  is 
highly  peaked  in  the  current  sheet  with  V2i  —  0  elsewhere.  Noting  the  relatively  flat  nfa 3) 
profile  across  the  current  sheet  (Figure  2,  third  panel),  we  see  that  it  is  this  sheared  V2i(x3) 
profile  that  provides  the  current  sheet  current.  This  is  in  contrast  to  the  Harris  model  in 
which  V2(x3)  is  uniform  with  a  peaked  ni(x3)  but  is  similar  to  the  symmetric  case  in  the 
Vo/vth  ^  1  regime  previously  studied  [Holland  and  Chen,  1993].  The  flow  velocity  for 
the  electrons,  V2e,  is  in  the  opposite  direction  as  determined  by  equation  (23).  The  third 
panel  shows  V3i,  which  is  perpendicular  to  the  current  sheet.  Recall  from  the  discussions 
in  section  3.4  that  Vu  —  VleB s  well  as  V&  =  V3e. 

An  interesting  effect  is  that  the  degree  of  asymmetry  in  the  asymptotic  distribution 
functions  alone  can  affect  the  current  sheet  thickness  and  the  total  current.  This  is 
a  nonlocal  and  purely  kinetic  effect  and  is  similar  to  the  process  (“kinetic  thinning”) 
whereby  the  current  can  be  enhanced  by  increasing  vn/vth,  or  by  increasing  the  field 
aligned  component  of  the  flow  velocity  of  the  particle  distributions  in  remote  source  regions 
[Harold  and  Chen,  1996]. 
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We  now  consider  the  forces  acting  on  the  equilibrium  current  sheet.  Recall  that  by 
construction,  the  electron  quantities  satisfy  the  electron  momentum  equation  (15).  For 
the  ions,  all  the  moments  are  obtained  as  averages  over  individual  ion  orbits,  and  they 
must  satisfy  equation  (14)  in  all  directions.  Thus,  there  must  be  nonzero  dQzki/dx^  for 
A;  =  1,  2,  and  3. 

In  order  to  verify  that  the  simulation  moments  satisfy  this  equilibrium  force  balance 
required  by  Vlasov  equation,  we  evaluate  for  the  x*.  direction  P3ki,  ra*  71*1^*14*,  and  the 
Lorentz  force  contribution  defined  by 

PLk{xz)  =  f  *  n<(x)(V<  x  B)  •  xfcdx,  (27) 

C  J—IjC 

where  x  is  the  integration  variable  and  k  =  1,  2,  3.  To  test  for  force  balance  in  the  a;* 
direction,  we  sum  these  contributions  and  define  a  vector  quantity  Ck(x3)  by 

Ck{x$)  =  —  [Q3ki(x3)  +  PLkfa)]  >  (28) 

where  the  normalization  factor  is  defined  by 


Pn  = 


(  p}3i  +  77l*7liV 3V3  + 


This  is  the  asymptotic  total  (kinetic  +  magnetic)  pressure  in  the  x3  direction  evaluated 
at  x3  =  —Lc.  The  significance  of  C*(x 3)  is  that  if  the  calculated  current  sheet  is  a 
valid  solution  of  Vlasov  equation,  Ck(x3)  is  equal  to  the  integral  of  the  fcth  component 
of  equation  (14)  across  the  current  sheet.  That  is,  the  simulation  current  sheet  satisfies 
equilibrium  force  balance  only  if  Ck(x 3)  is  a  constant  function  of  x3  for  each  k. 

For  the  Xi  direction,  Figure  4  describes  the  ion  moments  directly  calculated  from  in¬ 
dividual  ion  trajectories.  All  the  quantities  are  normalized  to  P Jy.  The  top  panel  gives 
the  P3u  element  of  the  ion  pressure  tensor,  which  is  found  to  be  comparable  in  magni¬ 
tude  to  the  directed  momentum  flux,  771*71*  V'li,  shown  in  the  second  panel.  The  fact 
that  the  two  sure  comparable  is  a  direct  result  of  the  strongly  sheared  flow  Vu  since  this 
strong  shear  results  in  a  large  change  in  the  directed  momentum  flow  across  the  cur¬ 
rent  sheet  in  the  x3  direction  comparable  to  the  change  in  P3i.  The  third  panel  shows 
Pli(x3)  =  -(ft/c)  f*lc  ni(x)V2i{x)Bndx- 

The  fourth  panel  shows  C\(x3)  for  the  Xi  direction.  We  see  that  C\  is  nearly  constant, 
showing  that  forces  in  the  x\  direction  are  well  balanced  due  to  the  nonzero  X3-gradients 
in  P31i  and  77i*7i*V3*Vi*. 

For  the  force  balance  in  the  x2  direction,  Figure  5  displays  normalized  P32*,  77ij  71^*1^, 
Pl 2,  and  C2.  The  pressures  here  are  relatively  small  compared  with  those  in  the  x\  and 
x3  directions.  The  main  point  here  is  that  the  sum  of  the  normalized  for  the  x2  direction, 
C2(x3),  is  constant  to  a  high  degree  of  accuracy,  showing  that  the  velocity  moments  of 
ion  trajectories  are  in  equilibrium  in  the  x2  direction. 

Figure  6  shows  normalized  P^u  771*7^ V^V^,  Pl3,  and  C3.  Comparing  this  figure  to 
Figures  4  and  5,  we  see  that  the  pressure  due  to  the  flows,  771*71* V^Vs*,  plays  a  lesser  role 
in  the  force  balance  along  x3  than  in  the  Xi  and  x2  directions.  In  the  x3  direction,  the 
plasma  pressure  P33*  (top  panel)  essentially  balances  the  Lorentz  force  term  Plz  (third 
panel).  In  the  second  panel,  we  expanded  the  vertical  scale  to  show  the  slight  asymmetry 
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(solid  line).  The  dashed  line  is  the  quantity  rii(x3)V3i(x3)  (arbitrary  units,  scaled  to  fit  the 
plot),  which  is  essentially  constant.  This  means  that  the  continuity  equation  V  •  ( riiV* )  = 
0  is  satisfied.  The  force  balance  in  the  Xz  direction  is  demonstrated  by  the  essentially 
constant  Cz(xz)  shown  in  the  bottom  panel. 

Figures  4-6  demonstrate  that  continuity  and  force  balance  for  the  ions,  equation  (14), 
are  well  satisfied  in  all  directions  by  the  computed  current  sheet.  The  off-diagonal  elements 
and  their  nonzero  gradients  in  xz ,  i.e.,  dQzki/dx3  0,  are  essential  for  this  balance. 

5.  A  Solar  Wind  Current  Sheet 
5.1.  WIND  Spacecraft  Data 

The  above  example  shows  that  the  full  width  of  collisionless  current  sheets  is  expected  to 
scale  with  the  ion  Larmor  radius,  with  pa  ~  100  km  in  the  solar  wind  at  1  AU.  Advecting 
at  the  SW  speed  of  400  km/s,  such  current  sheets  are  expected  to  appear  as  structures  of 
the  order  of  1  second  in  duration.  To  compare  our  model  results  with  in  situ  SW  data, 
we  use  the  high  time-resolution  (0.044  sec  sampling  rate)  WIND  MFI  magnetic  field  data 
[. Lepping  et  al.,  1995]  and  3-second  WIND  3DP  ion  data  [Lin  et  al.,  1995].  We  do  not 
include  the  high  resolution  (3  sec)  electron  data  since  the  electron  data  are  generally  not 
sufficiently  accurate  for  short  time  scales  due  to  photoelectron  contributions,  which  can 
only  be  corrected  for  by  processing  the  data  over  longer  periods.  For  this  reason  we  only 
use  the  average  ratio  of  the  electron  to  ion  temperatures,  Te/Tu  over  the  entire  simulation 
box  as  input  to  the  simulation.  The  WIND  MFI  magnetic  field  data  have  estimated 
measurement  uncertainties  of  ±0.1  nT.  The  major  measurement  errors  in  the  WIND  3DP 
particle  data  are  due  to  truncation  errors  in  processing  the  digitized  data  and  are  typically 
±2  km  s-1. 

In  Figure  7,  we  show  the  WIND  MFI  data  for  a  SW  magnetic  discontinuity  centered, 
appearing  as  a  sharp  reduction  in  \B\,  at  16:28:28.5  UT  on  1996  January  16.  The  open 
diamonds  are  the  3-second  time-averaged  WIND  MFI  data,  showing  the  magnitude  of  B 
and  its  three  components  in  the  GSE  ( x-y-z )  coordinates.  The  magnetic  field  is  actually 
sampled  at  0.044-second  time  resolution.  This  is  shown  by  the  jagged  line  overplotted 
on  the  3-second  data  points.  The  spin  period  of  the  spacecraft  is  3  seconds  so  that  the 
higher-resolution  data  show  noticeable  spin  effects.  The  3-second  data  shown  here  have 
been  averaged  over  the  spin  period  to  eliminate  these  effects.  The  figure  shows  that 
the  magnetic  discontinuity  is  now  defined  by  one  data  point  but  that  the  overall  width 
and  depth  of  the  magnetic  field  reduction  are  not  significantly  modified  by  the  averaging 
process.  However,  the  0.044-second  data  possibly  show  a  slightly  (~20%)  narrower  width 
than  may  be  inferred  from  the  3-second  average.  Based  on  the  Bx  and  By  components, 
we  see  that  the  sharp  reduction  in  B  is  due  to  a  field  reversal,  i.e.,  a  current  sheet. 

Figure  8  shows  the  magnitude  and  vector  components  of  the  ion  flow  velocity  V<  as 
measured  by  the  WIND  3DP  instrument.  The  data  are  shown  in  the  GSE  coordinates. 
We  use  the  highest  time-resolution  particle  data  available,  which  is  3  seconds  for  this 
event.  The  data  indicate  that  there  is  a  significant  amount  of  shear  in  the  ion  flows  across 
the  current  sheet. 

Overall,  we  find  that  an  asymmetric  current  sheet  is  most  clearly  reflected  in  anisotropic 
and  nondiagonal  pressure  tensors  and  asymmetric  ion  momentum  flows.  Such  asymmetries 
are  consistent  with  the  recent  CLUSTER  observations  [e.g.,  Runov  et  al.,  2006].  The 
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current  density  J2  does  not  show  significant  spatial  asymmetry  in  the  parameter  regime 
studied.  Malova  et  al.’s  [2007]  model  based  on  approximate  ion  dynamics  also  shows 
only  slight  spatial  asymmetry  but  did  not  discuss  the  properties  of  the  pressure  tensor  or 
asymmetric  flows. 

5.2.  Model-Data  Comparison 

For  the  purpose  of  comparing  the  data  with  the  model  results,  we  transform  the  3-second 
magnetic  field  data  into  the  coordinate  system  in  which  the  X3  axis  is  normal  to  the  plane  of 
the  field  reversal,  which  is  calculated  using  the  minimum  variance  method  [Sonnerup  and 
Cahill ,  1967].  The  transformation  is  given  by  equation  (5).  The  specific  rotation  angles 
for  this  SW  structure  are  (0y,9z>,0xn)  =  (—15,40,60)  in  degrees.  The  unit  vector  normal 
to  the  current  sheet  in  the  transformed  coordinate  system  is  X3  =  0.408x— 0.663y+0.627z 
in  terms  of  the  GSE  unit  vectors.  The  SW  velocity  is  Vsw  =  (-550,23,-15)  km/sec 
at  the  center  of  the  magnetic  discontinuity.  The  same  rotation  is  used  to  transform  the 
average  velocity  component  14*.  The  density  and  temperature  are  not  affected  by  this 
rotation. 

From  the  transformed  data  we  determine  the  asymptotic  values  of  B±,  nf,  Vxf ,  and 
7*  and  adjust  the  input  ion  distribution  functions  to  obtain  the  solution  that  has  the 
best  overall  fit  to  these  boundary  values.  The  specific  input  parameters  to  adjust  are  Uf , 
n±,  and  v 

Figure  9  provides  a  direct  comparison  of  this  model  current  sheet  (solid  curves)  and  the 
3-second  WIND  MFI  magnetic  field  data  in  the  transformed  frame  (open  diamonds).  Here 
and  in  the  subsequent  plots,  observation  time  is  converted  to  £3  using  the  x  component  of 
Vsw  and  the  projection  of  the  x  coordinate  onto  the  normal  to  the  current  sheet.  Thus, 
the  3-second  time  interval  corresponds  to  the  spacing  between  successive  open  diamonds. 
The  four  panels  show  that  the  simulation  results  are  in  excellent  overall  agreement.  In 
particular,  the  B\  component  (second  panel,  replotted  from  Figure  2)  clearly  shows  a 
hyperbolic  tangent-like  field  and  closely  fits  the  observed  field  reversal.  The  third  and 
fourth  panels  show  that  B2  and  B3  are  both  approximately  constant,  with  the  values 
B2  fa  0  and  B$  ~  1  nT.  Thus,  the  model  assumptions  regarding  the  magnetic  geometry 
are  well  satisfied  for  this  event.  The  overall  agreement  between  the  model  and  data  is 
good. 

The  top  panel  indicates,  however,  the  calculate  magnitude  of  the  asymptotic  magnetic 
field  is  ~15%  smaller  than  the  observed  values.  This  is  due  to  the  fact  that  the  calculated 
B\  is  about  15%  less  than  the  observed  values.  This  implies  that  the  the  total  current  in 
the  x2  of  the  simulated  current  sheet  is  less  than  that  of  the  actual  SW  structure.  While 
the  observed  reversal  in  Bi  is  well  replicated  by  the  model,  the  thickness  of  the  reversal 
region  is  not  well  resolved.  It  is  possible  that  the  actual  thickness  is  greater  than  that 
of  the  model  field  reversal,  perhaps  by  a  factor  of  1.5  to  2.  We  did  not  find  a  better 
fit  by  varying  the  input  parameters  given  the  form  of  the  asymptotic  ion  distribution 
function.  In  the  physical  system,  a  slight  thickening  of  the  current  distribution  may  occur 
due  to  instabilities  that  can  cause  diffusion  and  broadening  of  the  current  sheet  until 
nonlinear  saturation  occurs  (e.g.  Kelvin-Helmholtz  [Miura  and  Pritchett,  1982;  Pritchett 
and  Coroniti,  1984]).  Our  model  does  not  include  such  effects. 

Figure  10  displays  the  magnitude  and  components  of  the  transformed  velocity  data. 
Here,  the  transformed  velocity  has  been  Vsw  subtracted  out  so  that  the  data  and  model 
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results  can  be  compared  directly  in  a  common  reference  frame.  In  this  figure,  we  have 
overplotted  the  model  results  (solid  curves)  shown  in  Figure  3  on  the  WIND  3DP  data 
(open  diamonds).  The  top  panel  shows  the  magnitude  of  the  velocity,  followed  by  the 
individual  components,  Vu,  Vu,  and  V^.  The  particle  data  exhibit  more  fluctuations  than 
the  magnetic  field  data,  but  the  overall  agreement  is  good  with  a  few  notable  exceptions. 

For  Vu  (second  panel),  the  jump  in  Vu  from  —25  km/sec  to  20  km/sec  across  the 
magnetic  field  reversal,  the  degree  and  sense  of  asymmetry,  and  the  near  constancy  of 
Vu  on  either  side  of  the  jump  are  all  in  good  agreement  between  the  model  and  data. 
The  calculated  (third  panel)  is  nearly  zero  everywhere  except  around  the  midplane 
where  it  has  a  pronounced  peak.  This  is  in  disagreement  with  data.  The  observed  Vzi 
component  is  not  zero  and  does  not  have  a  simple  peak.  This  discrepancy  is  greater  than 
the  typical  error  in  velocity  measurements  of  of  the  order  of  ±2.  The  data  points  exhibit 
asymmetry  across  the  current  sheet,  varying  from  an  average  of  <  Vu  >—  5  km/sec  for 
x3  <  0  to  <  Vu  >~  8  km/sec  for  x3  >  0.  This  difference  in  the  measurement  values  is 
consistent  with  the  error  bars  but  is  probably  too  systematic  to  be  random  fluctuation. 
The  overall  deviation  from  Vu  =  0  is  also  greater  than  the  error  bars  throughout.  The 
most  notable  difference  is  the  pronounced  peak  in  V<u,  which  coincides  with  the  dip  in 
the  magnitude  of  velocity  |Vj|,  which  arises  from  the  reversal  in  Vx<.  The  bottom  panel 
of  Figure  10  shows  a  comparison  of  the  calculated  V®  component  with  the  observed  l/jj. 
Overall  agreement  in  magnitude  is  again  good.  In  addition,  the  observed  Vu  is  roughly 
constant,  as  is  the  simulation  V3i-  The  V&  data  also  exhibits  fluctuations  of  approximately 
±4  km/sec,  as  in  the  other  components,  consistent  with  but  slightly  greater  than  the 
estimated  measurement  errors. 

In  Figure  11  the  ion  density  and  temperature  profiles  for  the  WIND  3DP  data  (open 
diamonds)  are  compared  with  the  simulation  ion  density  and  temperature.  The  solution  is 
able  to  approximately  (within  ~20%)  match  the  observed  nf  and  T±.  The  model  results 
for  both  quantities  and  the  data  show  similar  profiles.  In  particular,  the  computed  rij 
jumps  from  n*  2.5  cm-3  to  n*  «  3  cm-3  across  the  current  sheet.  The  observed  n<  is 
nearly  constant  at  2.5  cm-3  for  x3  <  0  and  increases  to  3  cm-3  across  the  current  sheet, 
thereafter  slowly  decreasing  to  about  2.7  cm-3.  The  simulated  ion  temperature  and  the 
data  exhibit  similar  asymmetry  across  the  current  sheet,  but  the  observed  temperature 
jump,  A Ti  ~1  eV,  is  smaller  than  in  the  model  results  of  AT*  ~3  eV.  The  sharpest  jump 
in  temperature  occurs  on  the  positive  x3  side  of  the  current  sheet  in  both  the  data  and 
simulation  result.  The  plasma  &  of  this  SW  structure  is  /?+  «  1.5  in  S~  and  ^  w  30  in 
the  current  sheet. 

The  calculated  velocity  components  show  more  discrepancies  from  the  observed  data 
than  do  other  quantities,  as  described  abpve.  The  overall  offset  from  zero  in  the  Vu 
component  (third  panel,  Figure  10)  may  indicate  the  presence  of  an  overall  drift  velocity 
of  about  6  km/s  in  the  observed  V2»-  Had  a  constant  E\  component  been  included  in  the 
simulation  such  that  Vi e  =  cEi/Bn  «  6  km/s,  the  solution  would  have  been  uniformly 
offset  by  6  km/s  without  significant  changes  to  the  equilibrium  structure.  (A  heuristic 
discussion  of  the  effects  of  electric  field  is  given  in  Appendix  A.)  This  would  also  correspond 
to  nonzero  charge  separation,  which  is  insignificant  for  Va/c  -C  1.  The  presence  of  a 
uniform  E\  0,  however,  would  not  produce  the  observed  systematic  asymmetry  in 
Vu-  We  also  note  that  the  calculated  density  rii  and  temperature  7J  profiles  (Figure  11) 
show  narrow  peaks  of  about  20%  centered  at  the  current  sheet,  while  the  data  exhibit 
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no  apparent  peaks.  For  both  the  peaks  in  n*  and  7}  and  V2 »,  the  widths  are  comparable 
to  the  3-second  resolution  of  the  instrument.  It  is  possible  that  the  underlying  particle 
feature  was  not  resolved  by  the  3-second  resolution  data.  Comparing  the  Va  data  in 
Figure  10  and  n*  data  in  Figure  11,  we  see  that  if  the  observed  and  V2*  truly  had 
no  peaks,  their  product  would  also  have  no  peak,  implying  that  there  would  be  no  ion 
current  sheet.  Equation  (16)  shows  that  the  peak  in  J2  on  the  scale  of  pa  then  would 
all  have  to  come  from  the  electron  current  J2e-  For  the  simulation  parameters,  equation 
(23),  with  calculated  (•  •  *)M  «  0.6,  shows  that  Jn  ~  J2e •  Thus,  it  is  unlikely  that  V2(  had 
no  peak  corresponding  to  the  current  sheet,  supporting  the  possibility  that  the  peak  was 
not  resolved.  Another  possibility  can  be  traced  to  the  fact  that  the  width  of  the  observed 
magnetic  discontinuity  may  be  about  a  factor  of  two  wider  than  the  computed  width.  If 
this  is  true,  it  implies  a  wider  current  peak  carrying  approximately  the  same  total  current, 
yielding  a  peak  in  the  measured  V2i  that  is  broader  and  about  a  factor  of  two  lower  in 
height  than  what  is  shown  in  Figure  10. 

We  have  analyzed  several  other  thin  SW  current  sheets.  The  model  assumptions  appear 
to  be  well  satisfied  by  the  example  discussed  here.  As  such  it  is  a  good  example  of 
the  subclass  of  SW  discontinuities  with  Bn/B  ~  0.1,  B2  ~  0,  and  S  ~  pa.  The  other 
examples  have  the  same  basic  characteristics  (e.g.,  8  ~  pa,  asymmetry,  peaked  J2  but 
relatively  constant  niy  nondiagonal  pressure  tenser),  but  the  agreement  is  not  as  “clean,” 
indicating  that  the  model  idealization  may  not  be  as  applicable.  Overall,  the  model 
describes  the  other  examples  with  similar  degrees  of  agreement.  We  have  also  found 
examples  of  magnetic  discontinuities  whose  thicknesses  are  up  to  five  times  thicker  in 
terms  of  the  ion  Larmor  radius  than  the  examples  mentioned  here.  They  may  constitute 
a  different  class  of  currents  sheets. 

6.  Discussion  and  Conclusions 

We  have  presented  a  model  of  thin  current  sheets  that  occur  at  magnetic  discontinuities 
in  a  collisionless  plasma.  The  structure  is  1-D  but  is  allowed  to  be  asymmetric  about 
the  midplane.  The  model  is  valid  if  the  gradient  scale  lengths  in  the  plane  of  the  current 
sheet  are  longer  than  pn,  the  Larmor  radius  in  the  normal  component  of  the  field,  Bn 
and  represents  a  class  of  self-consistent  equilibria  in  which  off-diagonal  tensor  elements 
provide  force  balance  with  £3  =  0.  A  hybrid  kinetic  method  is  used  to  numerically  obtain 
time-stationary  solutions  of  the  Vlasov-Maxwell  equations  for  specified  ion  distributions 
in  the  asymptotic  regions.  A  /c-like  distribution  was  used  as  the  incoming  distribution 
function  for  the  ions.  The  current  J(x3)  consists  of  contributions  from  both  kinetic  ions 
and  massless  fluid  electrons.  The  scalar  potential  is  neglected,  and  quasi-neutrality  is 
assumed.  The  results  are  self-consistent  in  that  Ampere’s  law  J  =  (c/47r)VxB  is  satisfied. 
Equilibrium  force-balance  and  ion  continuity  are  explicitly  verified.  In  the  solar  wind  at 
1  AU,  the  ion  and  electron  contributions  to  the  current  are  comparable.  All  solutions 
we  have  obtained  converge  to  thicknesses  of  a  few  ion  Larmor  radii  (pa)  regardless  of  the 
initial  choice  of  Ba 0  and  80,  showing  a  posteriori  that  a  kinetic  model  including  the  full 
ion  motion  is  required.  The  present  model  is  a  significant  extension  of  such  past  works  as 
that  of  Lemaire  and  Burlaga  [1976],  which  did  not  include  the  nonlinear  particle  dynamics 
with  Bn^0  and  that  of  Holland  and  Chen  [1993],  which  was  limited  to  the  symmetric 
case,  both  of  which  excluded  strongly  sheared  flows.  As  in  Holland  and  Chen  [1993], 
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the  calculations  of  Malova  et  al.  [2007]  and  Mingalev  et  al.  [2009]  also  neglect  electron 
dynamics. 

The  results  show  that  equilibrium  current  sheet  solutions  do  exist  that  are  (1)  asymmet¬ 
ric  in  the  density,  current  density,  flow  velocity,  and  temperature  and  (2)  have  strongly 
sheared  flows  across  the  current  sheet.  We  modeled  in  detail  one  SW  current  sheet,  for 
which  high-resolution  data  are  available  for  both  the  magnetic  field  and  ion.  A  comparison 
of  the  resulting  current  sheet  structure  with  high-resolution  WIND  MFI  magnetic  field 
and  3DP  particle  data  shows  good  quantitative  agreement,  providing  a  well-constrained 
test  of  the  model.  While  we  cannot  prove  that  the  observed  current  sheet  discussed  here 
is  actually  in  equilibrium,  we  have  shown  that  the  observed  properties  are  consistent  with 
being  in  equilibrium.  We  stress  that  it  is  crucial  to  have  high-resolution  data  for  both  mag¬ 
netic  field  and  particles  for  comparison  with  the  model  results.  This  is  because  different 
current  sheet  configurations  with  different  particle  properties  can  have  similar  magnetic 
field  reversals  that  can  be  approximated  by  the  Harris  hyperbolic  tangent  profile.  The 
comparison  with  the  particle  data  provides  critical  constraints  on  the  results  in  terms  of 
matching  the  flows,  density,  and  temperature. 

We  have  found  additional  examples  of  possible  axisymmetric  current  sheets  with  thick¬ 
ness  of  the  order  of  pa  in  the  WIND  data,  but  the  model  fit  was  not  as  good  as  the 
example  discussed  in  the  present  paper,  with  a  comparable  degree  of  agreement  with  the 
magnetic  field  data  but  with  worse  agreement  with  the  particle  data.  The  data  do  not 
have  sufficient  temporal  resolution  to  identify  much  thinner  current  sheets,  if  any.  We  did 
not  find  any  structure  identifiable  as  current  sheets  having  thickness  significantly  greater 
than  pa. 

The  conclusion  that  asymmetric  collisionless  current  sheets  are  non-Harris-like  with 
strong  shears  in  the  average  speeds  and  nondiagonal  and  anisotropic  pressure  tensor  is 
important  because  these  plasma  properties  can  strongly  influence  the  stability  of  such 
current  sheets  and  play  a  critical  role  in  the  onset  of  reconnection.  For  example,  pressure 
anisotropy  can  strongly  affect  the  collisionless  tearing  mode  [Chen  and  Palmadesso,  1984; 
Chen  and  Lee,  1985;  Burkhart  and  Chen,  1989],  modifying  the  growth  rate  by  more  than 
an  order  of  magnitude  in  comparison  with  that  for  the  Harris  equilibrium. 

Note  that  a  model  current  sheet  is  chosen  to  best  fit  the  boundary  values  of  magnetic 
field  and  plasma  parameters.  The  properties  of  the  resulting  solution,  specifically  the 
thickness  and  magnitude  of  the  field  reversal  as  well  as  the  degree  and  sense  of  asymmetry, 
are  self-consistently  determined  by  the  physics  of  the  collisionless  ion  dynamics  and  fluid 
electrons  intrinsic  to  the  model. 

For  the  SW  current  sheets  we  have  examined,  only  a  moderate  degree  of  asymmetry 
is  required.  We  have  also  obtained  solutions  using  source  distributions  /+  and  f~  with 
much  greater  degrees  of  asymmetry.  Such  solutions  tend  to  have  greater  shear  in  the 
Vu(x3)  component  and  correspondingly  greater  current  ^(za)  than  current  sheets  with 
less  asymmetry  and  smaller  flow  velocity.  Thus,  the  degree  of  asymmetry  between  the 
two  spatially  separated  sources  remotely  affects  the  structure  of  the  current  sheet.  This 
nonlocal  effect  is  analogous  to  the  “kinetic  thinning”  process  in  the  symmetric  case  where 
the  current  can  be  enhanced  by  increasing  Vu/vth,  or  by  increasing  the  field  aligned 
component  of  the  flow  velocity  of  the  incoming  particle  distribution,  purely  kinetic  effect 
[Harold  and  Chen,  1996]. 
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It  is  apparent  from  our  sample  WIND  data  that  magnetic  discontinuities  and  the  asso¬ 
ciated  current  sheets  in  the  SW  can  be  thin  and  that  any  model  that  attempts  to  study 
these  structures  must  include  the  ion  Larmor-radius  scale  physics.  For  both  applications 
to  the  SW  and  the  magnetosphere,  the  central  physics  is  the  collisionless  motion  of  ions  on 
the  spatial  and  temporal  scales  of  ion  gyro  motion,  with  no  essential  limitation  imposed 
by  ambient  plasma  properties  such  magnetic  field  strength,  density,  and  temperature. 
The  good  agreement  with  the  high-resolution  WIND  data  and  the  previous  agreement 
demonstrated  for  the  magnetotail  current  sheet  suggest  that  the  basic  Larmor-radius 
scale  physics  underlying  our  model  current  sheet  is  scalable  to  current  sheets  in  regimes 
such  as  the  solar  and  stellar  coronae  and  other  astrophysics!  systems. 

Finally,  we  comment  on  the  electron  contribution  to  the  total  current.  We  have  treated 
electrons  as  a  massless  fluid  under  time-stationary  conditions  and  found  that  the  electron 
and  ion  contributions  are  comparable  for  Te  ~  7*.  In  the  quiet-time  magnetotail  where 
Te  -C  Tj,  our  model  would  result  in  j^el  <  l^il  according  to  equation  (23).  In  this 
situation,  diamagnetic  current  is  dominant.  However,  under  substorm  conditions  where 
there  is  a  strong  nonuniform  and  time- varying  electric  field,  it  has  been  noted  in  numerical 
simulations  [e.g.,  Pritchett  and  Coroniti ,  1995;  Hesse  et  al.,  1996]  and  magnetotail  data 
[. Mitchell  et  al.,  1990;  Asano  et  al.,  2003]  that  the  electron  contribution  can  dominate  the 
ion  contribution.  The  key  difference  is  that  these  simulations  and  observations  pertain  to 
driven  conditions  where  significant  charge  separation  may  build  up  in  the  current  sheet, 
rather  than  a  time-stationary  equilibrium  scenario.  It  is  interesting  to  note  that  even 
when  the  electron  contribution  to  the  current  is  dominant,  the  thickness  of  the  current 
sheet  is  still  of  the  order  of  the  ion  Larmor  radius.  The  present  work  shows  that  the  kinetic 
ion  physics  plus  hybrid  electrons  yields  thin  asymmetric  current  sheets  on  the  scale  of  ion 
Larmor  radius  in  agreement  with  observed  SW  structures  in  terms  of  both  magnetic  field 
and  particle  properties. 

Appendix  A 

In  this  appendix,  we  provide  a  heuristic  evaluation  of  the  deviation  (Sn/n)  from  exact 
charge  neutrality  due  to  the  asymmetry  in  the  pressure.  Equation  (8)  shows  that  Sn  can 
be  approximated  by 

_ 1  dEz  ^  El 

47re  dxz  4tt  epa  ’ 

where  the  definitions  6n  =  n,  —  ne  and  pa  =  vth/£la,  and  f20  =  eBa/miC  have  been  used. 
Here,  El  is  the  maximum  value  of  Ez,  which  occurs  at  xz  ^  6/2,  and  Ez  is  taken  to 
vary  over  the  thickness  of  5  ~  pa.  The  Ez  component  produces  an  E  x  B  drift  in  the  x% 
direction  in  the  frame  of  the  current  sheet  with  a  characteristic  value  of 


where  Ba  is  the  asymptotic  value  of  Bi  as  defined  in  section  1.  Using  the  above  equations, 
we  obtain 

*Ji„(Y±\2YLt 

n  \  c  J  vth 


(Al) 


(A2) 
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Since  Va/c  <C  1  for  solar  wind  conditions  and  we  consider  the  regime  &  vth  (analogous 
to  Vo/vth  &  1>  sections  1  and  2),  we  have 


n 


(A3) 


Charge  separation  is  second  order  in  Va/c  1  and  is  small  so  that  it  does  not  significantly 
influence  the  current  density  profile. 

A  more  precise  calculation  of  Sn  can  be  carried  out  by  directly  evaluating  £3.  Rrom  the 
ion  momentum  equation  (14),  we  find 


E  •  B  =  —(V  •  Qi)  •  B. 
en 

This  yields 


£3  = 


enBn 


dQzu  0  ,  dQsa  D 

01  +  — 1>3 


dx3 


dx3 


where  we  recall  E\  =  0.  From  equation  (28),  we  see 


dQzki  _  dPuc 
dx3  dx3  ’ 


(A4) 


(A5) 


(A6) 


where  k  =  1,  2,  4.  This  is  a  general  relationship  provided  the  definition  (27)  is  modified 
to  include  a  nonzero  E3.  For  the  special  subclass  of  equilibria  in  which  E3  =  0,  equations 
(A5)  and  (18)  yield 


dQ&i 

dx3 


1  8QsuBi  =  _enVxBv 


Bn  dx3 


(A7) 


The  constancy  of  C*  seen  in  Figures  4-6  demonstrate  that  force  balance  indicated  by 
these  relationships  is  well  satisfied  in  the  numerical  results. 

We  have  not  directly  evaluated  the  right  hand  side  of  equation  (A5)  because  of  the 
numerical  noise  in  the  derivatives  of  the  elements  of  Qi,  but  a  number  of  general  properties 
can  be  rigorously  inferred.  First,  we  note  that  £3  =  0  near  the  midplane  and  has  two 
extrema  away  from  the  midplane,  one  negative  and  one  positive  at  x3  ~  ±<5/2,  respectively. 
This  can  be  seen  as  follows.  For  the  first  term  in  the  square  brackets  in  equation  (A5), 
we  have  Bj  ~  0  at  S3  ~  0  because  of  the  field  reversal.  The  quantity  (dQ3u/dx3)  is 
peaked  in  the  current  sheet  because  it  is  the  flux  of  X\  momentum  in  the  x3  direction,  but 
because  Bi  is  essentially  an  odd  function,  the  product  (dQ3u/dx3)Bi  has  two  extrema 
of  comparable  magnitude  and  opposite  signs.  The  signs  of  the  extrema  depend  on  the 
asymptotic  distribution  functions.  In  our  example,  the  input  particles  predominantly 
enter  the  simulation  region  from  S~  so  that  the  extrema  are  negative  (positive)  for  x3  <  0 
(rr3  >  0).  This  can  be  seen  from  the  top  two  panels  of  Figure  4.  Second,  Q33i{x3)  is 
peaked  near  x3  =  0  because  this  is  the  sum  of  the  plasma  and  kinetic  energy  density 
(equation  (12))  so  that  dQ^/dxs  ~  0  near  the  midplane  with  two  extrema  of  nearly 
equal  magnitude  and  opposite  signs.  This  is  generally  true  and  is  also  the  case  for  our 
example,  as  can  be  seen  from  the  top  two  panels  of  Figure  6.  Thus,  £3  ~  0  near  the 
midplane  where  the  current  sheet  magnetic  field  is  minimum. 

An  important  implication  of  the  property  that  E3  ~  0  near  x3  =  0  is  that  no  significant 
particle  energization  due  to  £3  is  expected  near  the  midplane.  Furthermore,  the  fact  that 
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the  extrema  have  opposite  signs  means  that  they  tend  to  cancel  when  £3  is  integrated 
over  the  current  sheet  so  that  the  net  contribution  of  £3  to  the  total  current  is  small 
unless  the  degree  of  asymmetry  is  sufficiently  large. 

Equation  (A5)  shows  that  in  a  Harris-like  equilibrium  where  the  pressure  tensor  is  diag¬ 
onal,  Qzki  =  0  for  k  7^  3,  and  an  £3  electric  field  component  is  necessary  to  establish  force 
balance.  In  contrast,  in  the  class  of  equilibrium  discussed  in  the  present  paper,  the  forces 
in  the  X3  direction  are  balanced  by  the  off-diagonal  tensor  elements.  This  corresponds  to  a 
specific  set  of  boundary  conditions.  It  is  possible  that  different  equilibrium  configurations 
exist  corresponding  to  different  £3  components  determined  by  the  boundary  conditions. 
It  is  expected  that  the  inclusion  of  £3  ^  0  reduces  the  contribution  from  the  Q311  term. 
In  conjunction  with  equation  (A7),  this  implies  that  in  general 

|£3|  £  —  £1.  (AS) 

c 

This  provides  an  estimate  of  the  maximum  value  of  £3,  obtained  by  setting  the  other 
contribution  to  zero.  Thus,  the  E  x  B  drift  that  results  from  this  field  is  expected  to  be 
£3/ B\  ^  Vnj c  1. 

Note  is  that  the  ExB  drift  due  to  £3  is  the  same  for  both  the  ions  and  electrons  and 
has  no  contribution  to  the  total  current.  This  is  reflected  in  the  fact  that  equation  (22) 
contains  no  contribution  from  this  drift. 

In  Figure  10  (third  panel),  the  WIND  3DP  data  does  not  show  a  prominent  peak  in 
Vi i  near  the  midplane  of  the  current  sheet  while  the  model  result  does.  The  agreement 
between  the  data  and  the  model  can  be  improved  without  significantly  modifying  the 
structure  of  the  current  sheet  if  a  constant  drift  velocity  Vie  of  6  km/sec  is  added  to  the 
simulation  results.  This  would  imply  a  nonzero  E\  field,  corresponding  to  a  small  amount 
of  charge  separation  of  order  Vie/c  1. 

In  a  particle-in-cell  simulation  of  magnetic  reconnection  in  an  asymmetric  current  sheet, 
Pritchett  [2008]  found  a  strong  nonzero  electric  field  in  the  X3  direction  peaked  at  the 
current  sheet.  This  differs  from  the  equilibrium  £3  discussed  above.  If  such  a  peaked 
electric  field  existed  in  equilibrium,  one  would  expect  two  oppositely  directed  peaks  in 
Vu  across  the  midplane  of  the  current  sheet,  one  on  each  side.  The  transformed  WIND 
3DP  data  for  Vu  in  Figure  10  does  not  indicate  the  presence  of  such  features.  Comparing 
with  the  simulation  of  Pritchett  [2008],  this  indicates  that  the  observed  structure  was  not 
undergoing  reconnection  at  the  time  of  observation. 

There  is  another  ExB  drift  in  the  Xi  direction  that  is  produced  by  the  £1  component. 
This  is  distinct  from  the  above  drift.  In  particular,  £1  is  constant  in  a  1-D  system  (section 
3.2)  so  that  this  drift  is  constant.  In  Figure  10  (third  panel),  the  agreement  between  the 
data  and  the  model  can  be  improved  without  significantly  modifying  the  structure  of 
the  current  sheet  if  a  constant  drift  velocity  Vie  of  6  km/sec  is  added  to  the  simulation 
results.  This  discrepancy  may  be  a  manifestation  of  the  effect  of  a  nonzero  E\  field,  which 
is  neglected  in  the  present  calculation. 

Note  that  the  simulations  of  Hesse  et  al.  [1996J  and  Pritchett  [2008]  show  the  presence 
of  significant  electric  fields  only  during  reconnection.  The  assumption  of  negligible  scalar 
potential  for  equilibrium  considerations  needs  to  be  fully  validated  in  future  work  but  is 
consistent  with  these  simulation  results. 
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Finally,  we  note  that  the  constant  drift  velocity  in  the  Harris  model,  denoted  by  Vh, 
is  due  to  the  plasma  diamagnetic  drift  in  the  x 2  direction.  In  contrast,  the  drift  velocity 
discussed  here  is  produced  by  an  electric  field.  In  the  Harris  model,  S  =  {vth/V^Pai  so 
that  6  =  pa  when  %,  =  Vh-  In  contrast,  the  present  model  yields  S  =  pa  in  the  regime 
V-2E  <  vth,  consistent  with  the  SW  data. 
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Figure  1.  A  schematic  representation  of  a  current  sheet,  (a)  The  current  is  in  the 
x2  direction  (out  of  the  paper)  and  is  highly  peaked  between  the  dashed  lines  where  the 
“field-line”  curvature  is  maximum.  The  half-thickness  of  the  current  sheet  is  S.  The 
two  asymptotic  regions  (S+  and  S~)  far  away  from  the  current  sheet  (|x3|  S)  and 
the  associated  particle  distribution  functions,  f+  and  /“,  are  shown,  (b)  Schematic  of  a 
current  sheet  (box)  in  the  solar  wind,  x,  y,  and  z  are  aligned  with  the  GSE  coordinates, 
with  +x  pointing  from  the  Earth  toward  the  Sun,  +y  out  of  the  page,  and  +z  northward. 
The  X1-X2-X3  frame  is  co-moving  with  the  current  sheet,  with  £3  perpendicular  to  the 
plane  of  the  current  sheet.  The  asterisks  (*)  represent  the  entry  point  (right  *,  bottom 
surface)  and  exit  point  (*  on  left,  top  surface)  as  the  current  sheet  moves  past  an  observer. 


Figure  2.  Structure  of  a  self-consistent  model  current  sheet:  Bi  (x3)  is  the  magnetic 
field  component  in  the  £x  direction;  J2(x3)  is  the  current  density  in  the  y2  direction;  n<(x 3) 
is  the  ion  density;  T*(x3)  is  the  ion  temperature. 


Figure  3.  Components  of  the  average  ion  flow  velocity,  Vi*(x3),  V2*(x3),  V3*(x3),  in  the 
current  sheet  frame. 


Figure  4.  Terms  in  the  Xi  component  of  the  ion  force  equation  (14)  (E  =  0):  the 
pressure,  P3 M,  the  momentum  flux,  m*n*V3*Vii,  and  the  Lorentz  force,  Pn  (equation 
(27)),  all  normalized  to  Pn-  C\(xz)  is  the  sum  of  all  contributions.  Constant  C\  means 
that  force  balance  is  satisfied. 


Figure  5.  Terms  in  the  x2  component  of  the  ion  force  equation  (14)  (E  =  0):  the 
pressure,  P32i,  the  momentum  flux,  771*71* V^V^,  and  the  Lorentz  force  term,  P^  (equation 
(27)),  all  normalized  to  Ps-  C2(x 3)  is  the  sum  of  all  contributions.  Constant  C2  implies 
that  force  balance  is  satisfied. 
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Figure  6.  Terms  in  the  x3  component  of  the  ion  force  equation  (14)  (E  =  0):  the 
pressure,  the  momentum  flux,  rniniV^Vzi  (solid),  riiVu  (dashed,  arbitrary  scale  to  fit 

the  scale),  and  the  Lorentz  force  term,  Plz  (equation  (27)),  all  normalized  to  P^.  C3(x3) 
is  the  sum  of  all  contributions.  Constant  C3  implies  that  force  balance  is  satisfied. 


Figure  7.  WIND  MFI  data  of  the  magnetic  field,  B,  in  the  GSE  coordinates,  plotted 
with  t  =  0  at  16:28:28.5  UT,  1996  January  16.  The  open  diamonds  show  3  second  time- 
averaged  data.  The  error  bars  (±0.1  nT)  are  approximately  the  size  of  the  diamonds.  The 
data  is  sampled  at  0.044  sec  time  resolution  is  shown  by  the  thin  line.  The  jaggedness 
arises  from  the  spacecraft  spin  with  a  period  of  3  seconds. 


Figure  8.  WIND  3DP  data  in  the  GSE  coordinates.  Vi  is  the  average  velocity  for  the 
hydrogen  ions.  The  time  axis  uses  t  =  0  at  16:28:28.5  UT  1996  January  16.  Representative 
error  bars  are  shown,  approximately  ±2  km/sec  in  each  component.  For  Vi  and  Vn,  the 
error  bars  are  comparable  to  the  size  of  the  diamonds.  The  time  resolution  is  3  seconds. 


Figure  9.  A  comparison  of  WIND  MFI  magnetic  field  data  (open  diamonds)  in  the 
co-moving  (X1-X2-X3)  frame  of  the  current  sheet  and  the  simulated  equilibrium  current 
sheet  (solid  curves).  The  x3  axis  is  the  minimum  variance  direction  for  the  data  points 
and  is  normal  to  the  simulated  current  sheet.  The  field  reversal  occurs  in  B\. 


Figure  10.  A  comparison  of  WIND  3DP  data  (open  diamonds)  for  the  average  flow 
velocity  V*  of  ions  and  the  simulation  results  (solid  curves)  in  the  co-moving  frame  of 
the  current  sheet,  as  in  Figure  9.  The  Vu  component  shows  a  significant  amount  of  flow 
shear.  Representative  error  bars  are  shown. 


Figure  11.  A  comparison  of  WIND  3DP  data  (open  diamonds)  for  the  ion  density 
and  temperature  and  the  model  results  (solid  curves)  in  the  co-moving  frame  of  the 
current  sheet,  as  in  Figure  9. 


